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1 Introduction 

The study of forced oscillations emanating from a limit cycle is a classical problem in the theory of 
bifurcation. Around 1950 the basic method to deal with this problem was developed by Malkin in [TT] 
and this study was continued by Loud in [TU]. The state of the art before the contributions of Malkin 
and Loud can be found in the book by Lefschetz [5] . To describe the general framework we start with an 
autonomous system 

x = f(x) 

having a closed orbit T associated to a periodic solution Xo(t) with period T > 0. Notice that T is not 
necessarily the minimal period. The perturbation considered is 

* = f( x ) +eg{t,x;e) 

where g is periodic in t and its period is precisely T. The beginning of Malkin's method is the construction 
of a T-periodic function M — M(0) depending upon Xo(t) and g(-, -;0). The zeros of M are intimately 
linked to the possible bifurcations to T-periodic solutions for e > 0. Assuming some non-degeneracy 
conditions on xo(t) one can prove that if 6** is a non-degenerate zero of M (M(9*) = 0, M'(0*) ^ 0) then 
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the perturbed system has a family of T-periodic solutions satisfying 

x s (t) = x (t + 8*) + O(e), asejO. 

It is also possible to analyze the case of a zero of higher multiplicity (M(0») = 0, M'(9*) = 
0, • • ■ ,M( fc - x )(0*) = 0, mW(S«) ^ 0) but this requires long computations, see e.g. |10| and [§]. More 
recently a topological approach has been taken in [4j. A bifurcation exists as soon as 0* is a zero where 
M changes sign. The next step after the existence of bifurcating branches is the study of the stability 
properties. This was already considered in [11] , [10) and [B]. Assuming that T is an exponential attractor 
it can be proved that the bifurcating periodic solution is asymptotically stable when M'(0*) > and 
unstable when M'{6„) < 0. If 0* is a zero of a higher multiplicity, then the implicit function approach 
taken in [TU] and [BJ does not allow to detect bifurcation of stable periodic solution on the basis of the 
sign of M( fc )(0„) and some further computations have to be done. See in particular equations (3.5) in 6J 
and (4.23) in |10j . The purpose of our paper is to obtain a topological version of this result for increasing 
or decreasing zeros when the derivative of M at can vanish. In particular, we are interested in an 
unified answer which does not depend on the multiplicity of 0*. We will get a positive answer in the 
case of analytic systems. For this class of systems we will use a variant of Lyapunov-Schmidt reduction 
that will allow us to prove that if M is not identically zero then the number of T-periodic solutions is 
finite. This is inspired by the results of Nakajima and Seifert in [T2] and R.A. Smith in [15]. Once we 
know that T-periodic solutions are isolated we can talk about their topological index. This is just a 
localized version of the topological degree and the connections of this index with the stability properties 
of the corresponding solutions have been discussed in [7] [5] [5J Q2] . The computation of the index is then 
obtained via a result in the line of those in [1] . 

The rest of the paper is organized in three sections. In Section [2] we present some preliminary results 
on the autonomous system. The main Theorem as well as an example illustrating its applicability can be 
found in Section [3] This section also shows how to prove the main result via topological degree. Finally 
Section [3] is devoted to the proofs of three Lemmas previously employed. 

2 The autonomous system 

In this section we present some elementary facts about the non-perturbed system. They will be needed 
later in order to state our main Theorem. Let us start with the autonomous system 



defined on an open subset fl of R n . The vector field / : fl — > R™ is real analytic. 

Assume that Xo(t) is a non-constant periodic solution of {TJ with period T > 0. The associated 
variational equation is 



x = f(x) 



(1) 



y = f'(x Q (t))y. 



(2) 
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This is a T-periodic equation having the solution Xo(t). The Floquet multipliers are labelled as n±, ...,/i n 
and counted according to their multiplicity. It will be assumed that they satisfy 



Ml = 1) I A*2 1 < 1, \Vn\ < 1- 



(3) 



This condition implies that the closed orbit V = {xo(t), t S [0, T]} is an attractor (see pQ). The region 
of attraction is an open neighborhood of T which will be denoted by A. C fi. 

In view of the condition on the Floquet multipliers we know that the space of T-periodic solutions of 
((21) has dimension one. The same property must hold for the adjoint system 



(4) 



z = -f'(x (t)y z . 

The next result will provide an orientation in the space of T-periodic solutions of Q . 

Lemma 1 There exists an unique T-periodic solution z$(t) of satisfying 

(xo(t), Zo(t)) = 1, for any t £ M.. 

Proof. It is based on Perron's lemma [14] (see also [2], Sec. Ill, §12). This result says that if y(t) and 
z(t) are arbitrary solutions of @ and (01 then 

(y(t), z(t)) = constant. 

We will prove that if z\ (t) is a non-trivial T-periodic solution of (UJ) then 

(i(0),«i(0))^0. 

Since the space of T-periodic solution has dimension one this will complete the proof. 
To prove (0 we find a n x n matrix S such that 

/ 1 ... \ 





(5) 



S'- 1 F(T)5 



V 



.4 



where Y(t) is the matrix solution of J5]) with Y(0) — In and det(A — I) ^ 0. From the definition of S 
we have that its first column Si is an eigenvector of Y(T) corresponding to the eigenvalue fj,x = 1. In 
particular S\ is parallel to x(0). Consider the matrix X = (S2\--\S n ) composed by the remaining columns 
of S. From the definition of S and A, 

F(T)S = EA 

Next we apply Perron's Lemma to the solutions Y(t)Si and Zi(t), 



This implies 



(zx(Q),Si) = (zi(0), yCTJSi), i = 2,..,n. 



zi(0)*E = zi(0)*F(T)S = zi(0)*EA 
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Hence z\(0)*H(I — A) = and so z*(0)£ = 0. Now we can conclude that {5]) holds, for otherwise we 
should have z*(0)S = which is impossible if Z\{t) is non-trivial. 

□ 

As a simple example we consider the planar system 

X = (l — |a;| 2 ) x + z|a;| 2 a;, x = xi + 1x2 G C. 

It has the periodic solution aro(t) = e l * whose orbit T = S 1 attracts A = C — {0}. The period is T = 2Ntt, 
where N > 1 is an integer arbitrarily chosen. The variational equation along Xo(t) is 

y = (-1 + 2i)y + {-I + i)e 2lt y 

and has the Floquet solutions 

Vl (t) = x (t) = ie« y 2 {t) = e + i). 

In consequence /ii = 1 and ^ 2 = e _2T . The computation of zo(t) follows from the proof of Lemma[T] We 
know that 

(y 1 (t) 7 z (t)) = 1, (y 2 (t), Zo(t)) = constant = 0. 

The periodicity of e 2t j/2(i) and Zo(t) implies that this last constant must vanish. From these equations 
one obtains that 

z (t) = + 

3 Main result and an example 

Let us consider the perturbed system 

x = f(x)+sg(t,x,e), (6) 

where g : R x £1 x [0, £*] 1— > R n is continuous and T-periodic in t. We also assume that for each i e R 
the function g(t,-,-) has partial derivatives up to the second order with respect to (x,e) and these 
derivatives are continuous as functions of the three variables (t,x, e). The most important assumption 
on the regularity of g will be the analyticity with respect to x. This means that for each cc» 6 there 
exists r > such that if \\x — sc*|| < r then for j = 1, n 

gj(t,x,e)= 9a,j(t,e)(x-x*) a , t G R, £G[0,e*]. 

Here a — (ai, ...,a n ) is a multi-index and we employ the notation for powers x a = x" 1 ■ ... ■ x" n . The 
coefficients g a .j are continuous and T-periodic in t and the convergence in the above series is uniform in 
t and e. As in the previous Section the vector field / is real analytic on Q and this is enough to guarantee 
that the solutions of © depend analytically upon initial conditions once e and t have been fixed (see [5]). 
Again xo(t) is a non-constant T-periodic solution of (fTJ) satisfying ([3]). We consider the function 

M{6)= f (g(t,x (t + 6),0),z (t + 9))dt, 
Jo 
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where zq is given by Lemma [T] This function is T-periodic and real analytic and so it will have a finite 
number of zeros in [0, T[ unless it is identically zero. 

Given 0* G [0,T[ a zero of M, M(6*) = 0, we say that index(M,0*) = 1 if M{6) ■ (9 - 6>*) > when 
9 7^ 6* is close to 6*». When the inequality is reversed we say that index(M, 0») = — 1. In any other case 
we say that index(Af, 6**) = 0. 

Theorem 1 In the previous setting assume that M is not identically zero and let U be a bounded and 
open set satisfying 

TcU CU cA 

(Recall that T is the closed orbit associated to Xo(t) and A is its region of attraction). Then there exists 
Eq > such that if < e < £o the system (0j has a finite number of T-periodic solutions passing through 
U. Moreover, if 9* is a zero of M with index(M, 6*) ^ then there exists a T-periodic solution x e (t) of 
(0|) with 

x e (t) - x (t + 0*) -> as e[0, 

uniformly in t € R. This solution is asymptotically stable if index(M, 9*) = 1 and unstable if 
index(MA) = -1- 

To illustrate the result we consider the planar system 

x = (l- \x\ 2 )x + i\x\ 2 x + e(a(t) + b(t)x + c(t)x), (7) 

where x S C and a, b. c : M — > C are continuous and 27r-periodic. The autonomous system (e = 0) 
was already analyzed in the previous section and we can now construct the function M for xo(t) = e 4 *, 
Zo(t) = (1 + i)e lt and T — 2ir. A direct computation leads to the formula 

r2-rr 

M{9) = Re[ / (a(t)+6(t)e 4( * +e )+c(i)e- 4 (* +e ))(l-i)e- l(t+9) dt] = 27rRe 
Jo 

where a m , b m and c m refer to the Fourier coefficients of a, b and c, namely 

fm = ^ fj f(t)e~* mt dt. 

In principle Theorem [1] would provide information on a bounded region U whose closure is contained in 
C — {0}. However the specific properties of ([7]) will allow us to deduce global results. To illustrate this 
we first claim that for < e < 1 any 27r-periodic solution x{t) will satisfy 

max||z(t)|| <P+ := [1 + ||a||oo + ||6||oo + ||c||oo] 1/2 . 

Indeed if £» is an instant when m := max||x(i)|| = ||a;(t*)|| then the derivative ^||a;(i)|| 2 = 2 (x(t),x(t)} 
must vanish at i* . From the equation ([7]) we deduce that 

IM**)H 4 = ll^^*)!! 2 + £ (o(*») + b(U)x(U) + c(U)x(U),x(U)^ . 

It is not restrictive to assume that m > 1 and and by dividing the latter equality by m 2 the claimed 
estimate follows. Next we observe that x = is a 27r-periodic solution for e = 0. The variational 



(aie" i9 + b a + c 2 e- 2ie ) (1 - i) 
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equation is y = y with Floquct multipliers H\ = 11% = e 27r . A standard perturbation result guarantees the 
existence of some p_ £ (0, 1) such that, for small e, there is a unique 27r-periodic solution z e (i) satisfying 
max ||z £ (i)|| < p_. Moreover this solution is unstable since all Floquet multipliers are greater than one. 
Now we apply Theorem [T] in the region 

U = {x e C : p- < \\x\\ < p + }. 

The function M can be expressed as a trigonometric polynomial of the type 

M(6) = /3 + a cos{0 + 0) + cos2(d + tp), 

with (3 = 27rRe[6 (l - 2nai(l - i) = ae _i *, 2ttc 2 (1 - i) = -fe~ 2%v . Now it is clear that M is not 
identically zero if and only if 

|ai| + |Re(So(l-t)]| + |$!| >0. 

In such a case ([7]) has a finite number of 27r-periodic solutions passing through U, say N. From the 
above discussions we conclude that also the number of 2-7r-periodic solutions on the whole plane is finite, 
namely N + 1. When the function M does not vanish we obtain a uniqueness result: z e is the unique 
27r-periodic solution. When M changes sign we obtain at least two additional 27r-periodic solutions, one 
asymptotically stable and one unstable. Summing up, we observe that in this example the function M 
gives conditions for the existence and stability that are rather sharp. Notice also that the function M can 
have zeros of the type M(9 ) = M'(6 ) = M"(6 ) = 0, M'"(6 ) ^ and they lead to an asymptotically 
stable solution. 

Before the proof of the Theorem we will state three lemmas that will be proved in the next section. 
Our first preliminary result goes back to [TTJ page 387] and [TU] . It shows that the zeros of the function 
M are relevant for the location of T-periodic solutions. We shall say that a solution x(t) passes through 
a set S C M 2 if x{t) £ S for some real t. 

Lemma 2 Assume that Ek I is a given sequence and let Xk(t) be a T-periodic solution of 0) with e = Ek 
and passing through U . Then it is possible to extract a subsequence {xk(t)} and a number 6* £ [0, T[ such 
that M(0*) = and 

Xk{t) — xo(t + 6*) — > as k — > oo 

uniformly in t gt. 

For the next statements it will be convenient to employ the Poincare map V £ associated to © . Denoting 
by x(t;(,e) the solution of © satisfying x(0) — £, we notice that for small e and ( £ U this solution is 
well defined in [0, T]. This is a consequence of the theorem on continuous dependence since U is compact 
and for e — the solutions starting at U C A are globally defined in the future. This observation allows 
us to define 

V s :U^M. n , C^x(T;C,e). 

This map is analytic and its fixed points are in a one-to-one correspondence with the T-periodic solutions 
starting at U. 
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Lemma 3 Assume that 8q £ R is an isolated zero of M , then there exist e$ > and R > such that for 
any e £ (0, £o) the Poincare map V e of @) has at most a finite number of fixed points in Br(xq{6q)). 

The third preliminary result will establish a link between the index of the zeros of M and the fixed point 
index of the Poincare map. Results of this type were already obtained in [5j but we will present later an 
independent proof. The Brouwer degree of a map / on a domain f2 will be denoted by deg(/, Q). It is 
assumed that is open and bounded and / does not vanish on its boundary. 

Lemma 4 Assume that 8q is an isolated zero of M and V is an open neighborhood of xq(8q). Then there 
exist a number £* > and a family of open sets V e C R n , e £ (0, £*), satisfying 

X (6 )€V e , V E CV 

and such that 

deg(id — V e , V s ) = — index(0o, M), whenever e £ (0, £*). 
We are now in the position to prove Theorem [TJ 

Proof of Theorem [TJ If the function M does not vanish then © has no T-periodic solutions passing 
through U when e > is small enough. This is a consequence of Lemma [2] From now on we assume 
that M vanishes somewhere. Let T* > be the minimal period of Xo{t), so that T — kT* for some 
k = 1,2,... The function M has period T* and the sequence of zeros of M on [0, T*[ is denoted by 

o < e x < e 2 < ... < e m < t*. 

Another consequence of Lemma [2] is that for small e any T-periodic solution of ((6]) passing through U 

must remain close to the orbit T for all time. In particular we can assume that all T-periodic solutions 

passing though U have an initial condition corresponding to a fixed point of V e . 

Step 1. There exists e\ > such that if e £ (0, E\) then V e has a finite number of fixed points. 

Once again we apply Lemma [2] and restrict e so that all the fixed points are contained in some of the 

balls Bji(xo(6i)), i = 1, ...,m, where R is given by Lemma [31 The union of these balls contains all the 

fixed points of V e and we know by Lemma G3 that they contain a finite number of fixed points. 

We can also assume that R has been chosen so that these balls are pairwise disjoint. This will be 
employed later and it is possible since T* is the minimal period and so the points xo(9i) and xo(9j) are 
different whenever i =/= j. 

After this step we can define the index of a T-periodic solution passing through U. Assume that x(t) 
is such a solution for some e £ (0, £i). We can find an open set W C U such that x(0) £ W is the only 
fixed point of V £ lying on W. The index of x(t) is defined as 

yr(x) =deg(id-P e ,W). 

In principle this index could take any integer value but the condition ([3|) implies that 

7 t(z) S {-1,0,1}. (8) 
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This fact was already noticed by Krasnoselskii in [7] . We refer to [7] or [13] for the proof. 

Step 2. If x(t) is a T-periodic solution of (J5]) passing through U, then x(t) is asymptotically stable if 

jt(x) = 1 and unstable if 7t(^) 1- 

The condition and the continuity of the Floquet multipliers with respect to parameters imply 
the existence of a positive number a > such that if B(t) is a T-periodic and continuous matrix with 
1 1 -6 (i) 1 1 < a for all t then the system 

y = (f'(x (t))+B(t))y 

has Floquet multipliers fi*, ■ ■ ■ with [i\ positive and dominant and < 1 for i — 2, • • • , n. After a 
time translation we conclude that the same property holds for the more general class of systems 

y = (f'(x (t + 6)) + B(t))y, max||5(i)|| < a, B{t + T)=B{t). 

For small s any T-periodic solution passing through U has a variational equation in this class and so the 
Floquet multipliers have the structure described above. The conclusion of Step 2 is a consequence of [5] 
and [13]. 

Step 3. Assume that index(M, 8i) ^ 0. Then for any e £ (0, £2] the equation ^ has a T-periodic solution 
x with 

x(0) 6 B R (x o (0i)) and -f T {x) = -index(Af, t ). 
This is a consequence of Lemma [H Indeed we can find an open set V e C Bn(xo(9i)) with 

deg(id - V £ , V e ) = -index(Af, 6*,) 

and the additivity of the degree implies that 

m 

deg(id - V e ,V E ) = 7t(^), 

where the T-periodic solutions of ([6|) with Xj(0) £ 14. The conclusion follows from ([8|). 

Notice that the convergence of this periodic solution to Xq (t + 9i ) as e — >• is a consequence of Lemma [5] 
since the balls Bfi(xo(9i)) are pairwise disjoint. This completes the proof of the Theorem. 

4 Proofs of the Lemmas 

Proof of Lemma [2j We present a proof for completeness. Since Xk passes through U one can find 
Tfc G [0, T] such that Xk(Tk) G U. After extracting subsequences we can assume that 

T k ^T and a; fc (r fe )^C. 

Let x(t) denotes the solution of |T]) with initial condition x(t) = £. Since £ is a point in the region of 
attraction A we know that x(t) is well defined in [r, 00 [. By continuous dependence we know that Xh(t) 
converges to x(t) and the convergence is uniform on every compact interval where x(t) is well defined. 



In particular this applies to [t, t + T] and so x(r) = lima^ (r) = lima;^ (r + T) = x(t + T) . This implies 
that x(t) is a periodic solution of |T]). Since .4 is invariant for |T]) and x{t) G we deduce that the 
closed orbit associated to x must be contained in A. This implies that this orbit is precisely T and so 
there exists 9* £ [0, T[ such that x(t) = Xo(t + 0*). In particular Xk(0) — > xq(9*). It remains to prove that 
M(0*) = 0. To this end we consider the map 

< f(C,e)=^(C)-C, CeU, ee[0,e ]. 

This is a C 1 map and the derivative D^((, e) is an n x (n + 1) matrix. We claim that the rank of 
D$(xo(9*), 0) is strictly less then n. Otherwise the equation <&(£,£•) = should describe a curve in a 
small neighborhood of (xq(6*),0). However the set $ = contains the curve (a;o(6'),0) and also the set 
of points (xfe(0),£fe) accumulating on (xq(6*),0). Once we know that rank_D<i>(xo(0*), 0) < n, it remains 
to prove that 

rankD$(x (6l), 0) = n if M{6) =f 0. 
The partial derivative with respect to £ is the n x n matrix 

9 f $Oo(0), 0) = Y(T + 9)Y{6)- 1 - I n , 

where Y(t) is the matrix solution of Q with F(0) = /„. Again, the Fredholm alternative for linear 
endomorphisms is applied to deduce that 

Imd c $(x (9), 0) = [Ker {[Y(9)T 1 Y(9 + T)* - C ' " 

The kernel in the above formula corresponds to the initial conditions at time t = 9 of the T-periodic 
solutions of |4|. Hence it is spanned by z (9) and so 

hadc*(x (0),Q) = {v e M" : r)±zo(6)} . 

By differentiability with respect to parameters, the function y(t) = d e x(t,(,e) with £ = xq(9), e = 
solves 

y = f'(x (t + 6))y + g(t, x (t + 9), 0), y(0) = 0. 
A direct computation shows that 

- (y(t),z (t + 9)) = (g(t, x Q (t + 9), 0),z (t + 9)) 

and, integrating over the period, 

(y(T),z o (9))=M(0). 
When M{9) ^ the vector y(T) is not in the range of d^^(xo(9), 0) and so 

rank (d c $(x ((9), 0) |<9 £ <I>(> ((9), 0) ) = (n - 1) + 1 = n. 

□ 
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Proof of Lemma [3l It is based on a variant of the Lyapunov-Schmidt reduction. We divide it in four 
steps. 

1. The change of variables. The dominant eigenvalue of L — (To)' (xq(9o)) is /ii = 1 with eigenvector 
xq(9q). This eigenvalue is simple and so we can find a linear projection tt in W 1 satisfying 

7T 2 = 7T, nL = Ln, Kcr n = {Xx(9q); X G M}. 

This is so-called spectral projection and the hyperplane Y = Im(«rf — L) is invariant under L. 
Moreover, 

<r(L Y ) = {H2, -^Mn}, (9) 

where Ly : Y — > Y~ is the restriction of L to Y\ In the rest of the proof v denotes a generic vector lying 
in Y. 

Consider the map 

$ : (9, w) G M x Y" i — > x (9)+v G W 1 . 
This is an analytic function with partial derivatives at (9q, 0), 

d e n9 a ,0) = X (9 ), dY$(9 ,v) = idy. 

The Inverse Function Theorem implies that $ is a local diffeomorphism mapping (^OjO) onto xq(9q). In 
a neighborhood of this point we reduce the search of fixed points of V e to the equation V e o $ = $. More 
precisely we consider the equation 

V e (x (9)+v) = x {9)+v, \9-9 \<A, \\v\\ < A, (10) 

for some small A > 0. Notice that $ is independent of e and so A is uniform in e ^ 0. 

2. The auxiliary equation. The equation (fTU|) can be interpreted as a system in the unknowns 9 and v. 
As usual we apply tt and solve in v. This means that we look at the implicit function problem 

F{9 : v;e) := TrV E (x (9) + v) - nx (9) - v = 0. 

This function maps \9 — 9q\ < A, \\v\\ < A, e G [0, £*] into F and satisfies 

a,,F(0o,O;O) = L y - idy. 

From the condition (JH]) we deduce that the implicit Function Theorem is applicable and so we find r > 

and a : [9q — r, 9q + r] x [0, r] — > Y - such that 

7rP £ (x ((9) + a(6», e)) = 7rx ((9) + a(6», e). 

Moreover this is the only solution of F(9,v;e) = in some ball ||i>|| < R. The function a is of class C 1 
and analytic with respect to 9. Due to the uniqueness of a we have a(9, 0) = for any 9 G [6o — r, 9q + r] , 
which can be combined with the smoothness of a to find a number fi > such that 

\\a(9,e)\\^efi for any 9 G [9 - r, 9 + r] , ee[0,r]. (11) 
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In this process it can be necessary to reduce the size of r. 

3. The bifurcation equation. Assume that x(t;S, e) is a T-periodic solution of © with S close to xo(#o) 
and e small and positive. We know from the previous steps that the initial condition can be expressed as 

S = x (Q) + a(B,e) 

for some G £ [9q — r, 9q + r] . Our next task is to show that O must be a zero of the function 

M e {6) f (b e (t,6),z (t + 8))dt 
Jo 

with 

b e (t, 9) := g(t, x(t, £, e), e) - - £ [f{x{t, £, e)) - f(x (t + 9)) - f(x (t + 9)) ■ (x(t, £ e) - z (t + 0))] 
and 

C-a;o(0) + a (0,e). 

By construction y(t) = x(t, 3, e) — xo(t + O) has to be a T-periodic solution of the linear equation 

y = f'(x (t + Q))y + eb s (t,e). 
The Fredholm alternative implies that O is a zero of M e . 

4. Conclusion: the role of analyticity. In view of the previous steps it is enough to show that the function 
M e has a finite number of zeros in [8q — r, 8q + r] for small e. 

Since a(9,0) =0we obtain by continuous dependence that 

b £ (t,6) -> g(t,xn(t + 6),0) as e^O 

uniformly in t E [0,T] and 9 G [9q — r,9o + r]. Indeed we also need to use that / is smooth and the 
estimate (ITT|) . This is required to prove that the term related to / goes to zero. Also the differentiability 
with respect to initial conditions and parameters plays a role here. 

The function M e converges to M as £ — > uniformly in 9 e [9o — r, 9q + r] . We are assuming that M 
is not identically zero and so the same must happen to M e for small e. Since M £ is analytic we conclude 
that it has a finite numbers of zeros in [9q — r, 9o + r]. This is valid for e G]0, £o[ with eo > sufficiently 
small. 

□ 

Remark The standard Lyapunov-Schmidt reduction for the equation P s (£) = £ would start with the 
splitting 

£, = tix (9 ) + v, v eR, veY, 

and considering the system 

TTp £ (rixo(9o) +v) = v 
(id - ir)T> e {rix a (9 ) + v) = r)x (9 ). 
Instead of this we are considering a sort of nonlinear splitting induced by the change of variables of 
Step 1. The advantage is that our bifurcation equation leads directly to M{9) = as e { 0. The same 
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approach is taken by Hale and Taboas in [3], but they prefer to work in an infinite dimensional framework. 

Proof of Lemma [4j First we pick up any n — 1 linearly independent solutions yi, ...,y n _i of 
whose initial conditions at 9q satisfy (t/i(#o)> zq(9q)) = 0. Next we consider the n x (n — 1) matrix 
Yi(0) = (yi(6»)| . . . |y„_i(6»)) and notice that 

Y 1 (9 + T) = Y 1 (9)A e (12) 

where A# is a (n — 1) x (n — 1) matrix with eigenvalues /Uz, . . . , /X?,. To verify this it is enough to 
observe that the hyperplane Vg spanned by y%{9), . . . , y n -i(9) is invariant under the monodromy operator 
Mg : y(9) i— » y(# + T). This is a consequence of Perron's Lemma. The eigenvector of Mg associated to 
/ii = 1 is £o(0) and does not belong to Vg. In consequence the restriction of Mg to Vg has eigenvalues 
fj,2, . . . , fj, n . The matrix is precisely the representation of this restriction with respect to the basis 
yx{6), . . . , y n -\{9). This property of the matrix Y\(9) will be employed several times. First we will 
employ it to evaluate the topological degree of the auxiliary map 

$ e (6, C) = ~eM(9)x (9) + {Y 1 (9) - Y x {6 + T))( 

with respect to fia := (9o — 3,6o + x Bs(0), where <5 > is a small number and B$(0) is an open 
ball in M" _1 . We will impose several restrictions on the size of 8, the first being that M has no zeros on 
[9q — S, 8q + 5} other then 9q. Notice that by linear independence the equation Q £ (9, £) = in £1$ splits 
as M{6) = and (Yi(0) - Y 1 (9 + T))( = 0. Then 9 = 9 Q and from the identity $F2S) we deduce that 
C = 0. Thus the degree we want to compute is well defined and does not change if we replace Qs by any 
sub-domain containing (#o, 0). For the effective computation we diminish S > in such a way that $ £ is 
linearly homotopic to the vector field 

$(0,C) = -M(9)x (9 ) + (Y 1 (9 ) - Y 1 {6 () +T))C 

for e > sufficiently small so that dcg($ e ,fi 5 ) = deg($,rj 5 ). The matrix S = (x (9 )\Yi(9 ) -Yi(9 + T)) 
is non-singular and the map $ can be expressed as So [(— M) x id]. By the theorems on the evaluation of 
the topological index of a composition of vector fields (see e.g. [8], Theorem 7.1), of a product of vector 
fields (see e.g. jS], Theorem 7.4) and of a linear vector field (see jS], Theorem 6.1) we have that 

deg($, n s ) = index((6» , 0), $) = indcx(0, S) ■ indcx((6» , 0), (-M) x id) = -signdet S ■ index(6» , M). (13) 

Another restriction on 5 that will be useful later is related to the map ip(9,Q = Xq(9) + Yi(#)£- This 
map must be a diffeomorphism from tt$ onto its image and C V. Notice that this is possible since 

detip'(9 ,0) = det(x (9 )\Y 1 (9 )^0. 

Our next step is to show that the vector fields 

T £ (9,C) = (id-V e )(xo(9)+Y 1 (e)0 

and $ £ are homotopic on a sub-domain of Q$ for e > sufficiently small. Let x(t] 9, C, s) be the solution 
of (J6j) satisfying x(0) = xq(9) + Y\(9)(. The Taylor expansion leads to 

x(t,9,C,e) = x {t + 9) + Y 1 (t + 9)( + s f Y(t + 9)Y(s + 9)- 1 g(s 1 x {s + 0),O)ds + 0(e 2 + \\(f), 

Jo 
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where we recall that the matrix Y(t) was defined in Section O This expansion is obtained by computing 
the derivatives with respect to C and £ and applying the formula of variation of constants. The matrix 
Y*(t + 8)~ 1 Y*(9) is fundamental at t = 6 for the adjoint system and so 

z (t + 9) =Y*{t + 6)- 1 Y*(d)z Q (6). 

From the periodicity of z$ we deduce that 

Y*(T + 9)z (9)=Y*(9)z (9). 

Thus, 

^jf Y(T + e)Y- 1 (s + e)g(s,x (s + d),0)ds,z (9)\ = 

= [ (g(s,x (s + 9),0),z (s + 9))ds = M(9). 
Jo 

In consequence, 

(T e {9, C), zq(0)) = -sM{9) + 0(e 2 + ||C|| 2 ), T e {9, Q = (Y^ff) -Y 1 (6 + T))C + ej(9) + 0(e 2 + ||C|| 2 ), 

where 7 is defined by an integral. Perhaps after a new reduction of the size of 6 we can find a positive 
constant A such that 

max \((Y 1 (9)-Y 1 (9 + T))C,y k (9))\>A\\C\\, for every C G M"" 1 and \9 - 9 \ < 6. 

k= l,...,n— 1 

To justify this assertion we notice that, by continuity, it is enough to check it for 9 = 9q and in this 
case it follows from (TT^|) since (Yi(6>o) — ^1(^0 + T)) — Yi(9q)(I — Ag ) and (/ — Ag ) is non-singular. 
From now on the number S will be kept fixed. We are going to compute the degree of J- £ on the set 
W s = {{9,Q ■ \0-9 Q \ <5, I ICI I < £ 2/3 }- The boundary of W £ is composed by A x : 9 = 9 ±5, ||C|| < £ 2/3 
and A 2 : \9 - 9 \ < S,\\(\\ = e 2 / 3 . On A x we observe that for e small enough 

sign (F e (0, 0,z (9)) = -signM(^), with 9 = 9 ± 6. 

On A2 we claim that for some fc = 1, . . . , n — 1 (depending on £), 

sign {F e (9, C), y k {9)) = -sign {^(9) -Y 1 (9 + T))(, y k (9)) . 

Indeed, from the expansion of T e we find that for each k 

(F e (9,O,y k (0)) = - ((Yi(8) - Y 1 {9 + T))C,y k {9)) + 0{e). 

For some k, \ ((Yi(9) — Yi(9 + T))(, yk(6)) \ > Ae 2 / 3 and this term is dominant, leading to the coincidence 
of the signs. Summarizing, for (9, £) £ dW e the vectors $ e (9, C) and !F £ (9X) do not point in opposite 
directions and, therefore, the vector fields $ 6 and T e are linearly homotopic on W £ (see jH theorem 2.1]). 
By excision, 

deg(^ e , W e ) = deg($, Sis) - -signdet S ■ md(9 ,M). (14) 
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To finish the proof we define V £ = tp(W e ) and observe that (id — V e ) o i/j = T e on W e . The theorem on 
the degree of the composition implies that 

deg(td - V e , V e ) ■ deg(^ -x (e ),W £ ) = deg(^ e , W e ). 

For instance, Theorem 7.2, Formula 7.6 in [5] is applicable since dV £ = ip(dW e ), V e is connected and 
xa(0o) G V e . By the linearization theorem for topological degree (see e.g. [8], Theorem 6.3) we have that 

deg^ -x Q (0 Q ),W e ) = sign det 1>'(9 ,0) = sign det(i (^o) 1^1(^0)). (15) 

The conclusion of the Lemma follows from these last identities and (TT4"|) because 

signdet(x (9 )\Yi(6 )) = sign det S. (16) 

To prove this claim we consider the family of matrices 

*i(0o) - Ali(0o +T) = YtfoW - XA 9o ), X G [0, 1], 

where once again we have used (|12p. For A = and A = 1 we obtain the second blocks of the matrices 
appearing in the identity (I16| . The eigenvalues of Ag are /Z2, . . . , all of them with modulus less than 
one. Hence 

dst(±o(Ba)\Yi(0 o ) - \Yt(e + T)) + 

for all A G [0, 1] and so the sign of this determinant is independent of A. The identity (fTH|) expresses this 
fact for the extreme values of A. □ 
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